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ABSTRACT 

Absorption-line spectroscopy of 23 background QSOs and numerous background galaxies has let us 
measure the spatial distribution of metals and neutral hydrogen around 1044 UV-selected galaxies at 
redshifts 1.8 ;$ z £ 3.3. The typical galaxy is surrounded to radii r ~ 40 proper kpc by gas that has 
a large velocity spread (Aw > 260 km s _1 ) and produces very strong absorption lines (New 3> 10 14 
cm~ 2 ) in the spectra of background objects. These absorption lines are almost as strong as those 
produced by a typical galaxy's own interstellar gas. Absorption with an average column density of 
New — 10 14 cm~ 2 extends to ~ 80 kpc, a radius large enough to imply that most strong intergalactic 
CIV absorption is associated with star-forming galaxies like those in our sample. Our measurement 
of the galaxy/CIV spatial correlation function shows that even the weakest detectable CIV systems 
are found in the same regions as galaxies; we find that the cross-correlation length increases with CIV 
column density and is similar to the galaxy auto-correlation length (tq ~ 4/i _1 Mpc) for New ~ 10 12 5 
cm~ 2 . Distortions in the redshift-space galaxy-CIV correlation function on small scales may imply 
that some of the CIV systems have large peculiar velocities. Four of the five detected OVI absorption 
systems in our sample lie within 400 proper kpc of a known galaxy. Strong Lyman-a absorption is 
produced by the intergalactic gas within 1ft, -1 comoving Mpc of most galaxies, but for a significant 
minority (~ 1/3) the absorption is weak or absent. This is not observed in SPH simulations that omit 
the effects of "feedback" from galaxy formation. We were unable to identify any statistically significant 
differences in age, dust reddening, environment, or kinematics between galaxies with weak nearby HI 
absorption and the rest, although galaxies with weak absorption may have higher star-formation rates. 
Galaxies near intergalactic CIV systems appear to reside in relatively dense environments and to have 
distinctive spectral energy distributions that are characterized by blue colors and young ages. 
Subject headings: galaxies: formation — galaxies: high-redshift — intergalactic medium — quasars: 
absorption lines 



1. INTRODUCTION 

A tremendous amount of gravitational energy is re- 
leased when a massive galaxy forms, roughly 10 62 erg 
from its billion supernova explosions and another 10 62 
erg during the assembly of its central 10 8 -M Q black hole. 
Although 99% of the supernova energy is carried away by 
neutrinos, and a similar fraction of the black hole radia- 
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tion escapes the galaxy, the remaining energy is enough 
to unbind most of the galaxy's gas. The fate of the galaxy 
depends on what happens to the energy. If it is absorbed 
by sufficiently dense gas, it will be converted into radi- 
ation by two-body processes and will stream harmlessly 
away; but otherwise it will set in motion an enormous 
blast-wave analogous to those produced by nuclear ex- 
plosions on earth. These blast-waves, tearing through 
galaxies, stripping away material, plowing outwards, and 
coming finally to rest in distant intergalactic space, are 
believed to have significantly altered the evolution of the 
baryonic universe. 

Indirect evidence for their existence comes in two 
forms. First, if the gas near galaxies had not been heated 
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by some source, disk galaxies would be smaller than they 
are (Weil, Eke, & Efstathiou 1998) and the correlation 
between galaxy clusters' X-ray temperature and luminos- 
ity would have the wrong shape (Kaiser 1991; Ponman, 
Cannon, & Navarro 1999). Second, something quenches 
most galaxies' star formation at some point in their his- 
tory, and blast-waves seem the best candidate. Without 
them small galaxies would be too numerous (e.g., Cole 
et al. 1994), big galaxies would be too blue (Springel, Di 
Matteo, & Hernquist 2005) and too large a fraction of 
baryons would have been turned into stars (e.g., White 
& Rees 1978, Springel & Hernquist 2002). 

Direct evidence remains elusive, however. Although 
local star-forming galaxies are surrounded by outflowing 
gas (e.g., Heckman et al. 1990; Martin 2005), typically 
their observable outflows extend to only a few kpc, far 
smaller than their virial radii. It is unclear how far these 
flows will advance before they stall or how severely they 
will affect the evolution of the galaxies and their sur- 
roundings. The presence of metals in the intergalactic 
medium might be taken as proof that the blast-waves 
spill out of some galaxies' halos, but even that is con- 
troversial. Intergalactic metals could also have been pro- 
duced by the first generation of Population III stars or 
stripped from galaxies by other processes (e.g., Gnedin 
1998). 

Five years ago, our group began to search for direct 
evidence for large-scale outflows around galaxies at high 
redshift. High-rcdshift galaxies were attractive targets 
because they have larger star-formation rates and (possi- 
bly) shallower potentials than their nearby counterparts. 
In addition, the higher density of the universe in the past 
is a significant advantage, because it strengthens many 
of the absorption lines that can be used to detect inter- 
galactic gas. We identified three observations that would 
qualify, in our view, as direct evidence for large-scale out- 
flows: gas at large galactocentric radii moving outwards 
at greater than the escape velocity; a strong association 
of intergalactic metals and galaxies; and disturbances, 
near galaxies, to the lattice-work of intergalactic HI that 
pervades the high-rcdshift universe. 

We searched for this evidence by systematically map- 
ping the relative spatial distributions of star-forming 
galaxies and intergalactic gas. Early results were pre- 
sented in Adelberger et al. (2003). This paper is an 
update. After describing the current data in § [21 we 
review the status of our search. Section [21 shows that 
the outflowing material responsible for the blue-shifted 
"interstellar" absorption lines in the galaxies' spectra of- 
ten lies at a radius of ~ 20-40/i -1 proper kpc. Weaker 
absorption lines extend out to radii approaching 200/i~ 1 
kpc. Section0]describes the galaxy-CIV cross-correlation 
function and reinforces the idea that detectable inter- 
galactic metals tend to lie near actively star-forming 
galaxies. Section OD shows that most galaxies are sur- 
rounded by significant amounts of intergalactic HI but 
that galaxies with little nearby HI are significantly more 
common than would be expected in the absence of winds. 
SectionEldiscusses the characteristics of galaxies that ap- 
pear to be associated with unusually weak Ly-a or un- 
usually strong CIV absorption. As we discuss in § 
these observations may be suggestive of galactic winds, 
but they do not rule out other alternatives. Our conclu- 
sions are summarized in § [HI Throughout the paper we 



Table 1. Observed fields 



Field 


Af2 a 


Nlris 


Nnirspec" 


Q0000 


3.4 x 4.5 


18 


1 


Q0201 


6.8 x 7.4 


25 


3 


Q0256 


7.3 x 6.6 


47 


1 


Q0302 


6.5 x 6.9 


44 


1 


Q0933 


8.1 x 8.0 


64 





Q1305 


11.8 x 11.0 


80 


1 


Q1422 


7.2 x 14.2 


113 


1 


Q1623 


16.1 x 11.6 


213 


31 


Q1700 


11.5 x 11.0 


89 


16 


Q2233 


8.2 x 8.6 


16 


1 


Q2343 


22.5 x 8.5 


191 


19 


Q2346 


10.9 x 11.0 


50 


7 


SSA22a 


7.8 x 8.4 


61 


8 



a Area imaged in U n G1Z and observed spectroscopically, arcmin 2 
b Number of galaxies with LRIS redshift z > 1 
c Number of galaxies with NIRSPEC redshift z > 1 



assume a cosmology with Q^j — 0.3, f^A = 0.7, h = 0.7. 
2. DATA 

The galaxies in our analysis were taken from the red- 
shift surveys of Steidel et al. (2003; 2004). These surveys 
targeted galaxies with magnitude 1Z £ 25.5 whose U n G1Z 
colors suggested that they had redshift 2 ^ z £ 3 (Adel- 
berger et al. 2004; Steidel et al. 2003). We restricted 
our analysis to galaxies in the 13 survey fields that con- 
tained one or more background QSOs. Although these 
fields include those previously analyzed by Adelberger et 
al. (2003), much of the present analysis relies on near-IR 
nebular redshifts (see below) that are only available in 
other fields. As a result, the overlap is small between the 
samples analyzed here and in Adelberger et al. (2003). 

Our 13 fields are scattered around the sky and have 
typical area Ail ~ 100 arcmin 2 (Table GJ. Redshifts for 
a total of 1044 objects in these fields were measured from 
low-resolution (~ 10A) multislit spectra taken between 
1995 and 2004 with LRIS (Oke et al. 1995; Steidel et 
al. 2004) on the Keck I and II 10m telescopes. The 
typical spectroscopic exposure time was 1.5 hours per 
multislit mask, although a subset of the objects was ob- 
served for more than ten times longer (Shapley et al. 
2005, in preparation). The redshifts obtained from these 
spectra are imprecise, because they rely on features (Ly- 
a emission and various far-UV interstellar absorption 
lines) that are redshifted or blueshifted relative to the 
galaxy's stars. Erb et al. (2005, in preparation) used 
NIRSPEC (McLean et al. 2000) on the Keck II tele- 
scope to obtain more-precise redshifts from the nebular 
emission lines ([OII]A3727, Ha, H/3, [OIII]AA4959, 5007) 
of 90 galaxies in this sample, many near the sightlincs 
to background QSOs, and we adopt their redshifts wher- 
ever possible. If nebular redshifts were unavailable, we 
estimated the stellar redshifts with the relationships de- 
fined in the following paragraph. Precise redshifts were 
required for much of the analysis, however, and in these 
cases we ignored galaxies that had not been observed 
with NIRSPEC. 

We measured the relationship between redshifts from 
Lyman-a emission (zLya), interstellar absorption (zism), 
and near-IR nebular emission (z ne b) for the full 138- 
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Fig. 1. — Velocity differences between the Lyman-a emis- 
sion line, interstellar absorption lines, and near-IR nebular lines 
for galaxies in the NIRSPEC sample. Each point represents one 
galaxy; shaded regions show the mean velocity difference ± stan- 
dard deviation of the mean for two redshift bins. The value of 
Spearman's rank correlation coefficient and its significance are 
r a = 0.16, P = 0.19 and r s = -0.197, P = 0.03 for the data 
in the left and right panels, respectively. The weak apparent cor- 
relations of velocity offset with redshift are therefore significant at 
only the 1-2<t level. 



object NIRSPEC sample of Erb et al. (2005, in prepa- 
ration), and used the results to help assign systemic red- 
shifts to galaxies when NIRSPEC redshifts were unavail- 
able. Figurenshows the observed velocity offsets for the 
galaxies in the full NIRSPEC sample. The least-squares 
first-order fits to the data can be written 

^ncb = ZLy* - 0.0033 - 0.0050(z Lytt - 2.7) (1) 
for galaxies with no detected interstellar absorption lines, 

^ncb = ^ism + 0.0022 + 0.0015(zi SM - 2.7) (2) 

for galaxies with no detected Lyman-a emission, and 

z„cb = z + 0.070Az - 0.0017 - 0.0010(z - 2.7), (3) 

with z = {z Lya + Zism)/2, Az = z Lya - z ISM , for galaxies 
with detectable Lyman-a emission and interstellar ab- 
sorption. The rms scatter around the three relationships 
is a z = 0.0027, 0.0033, 0.0024, respectively. In each equa- 
tion the final term accounts for the weak apparent trend 
of smaller velocity offsets at lower redshifts. Figure |2] 
shows the expected error distributions for the galaxies 
that lack near-IR nebular redshifts, calculated by com- 
bining the observed error distributions for equations^l2J 
and |3 after weighting by the number of galaxies in each 
class (233 with Lyman-a emission only, 425 with inter- 
stellar absorption only, 296 with both). For comparison, 
repeated observations of galaxies in the NIRSPEC sam- 
ple suggest that these galaxies' redshifts have a random 
uncertainty of 60 km s _1 . 

The criteria used to select galaxies for spectroscopy 
were irregular. We favored objects with 23 < 1Z < 24.5 
and objects near any known QSOs, and strongly disfa- 
vored objects far from the field centers. Near QSOs we 
often observed objects whose colors satisfied any of the 
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Fig. 2. — Estimated histogram of velocity errors Av = v es t — vtmc 
for galaxies without near-IR nebular redshifts. 



U n G7Z selection criteria of Steidel et al. (2003) or Adel- 
bergeretal. (2004). Far from QSOs our spectroscopic se- 
lection tended to be prejudiced towards one or the other 
color-selection criterion, "BX/BM" or "LBG," depend- 
ing on the field. This haphazard approach imposes a 
complicated selection bias on our data. Rather than try 
to correct it with elaborate simulations, we restrict our 
analysis to quantities that are not affected by irregular 
angular sampling or the assumed redshift selection func- 
tion. Further details of the galaxy observations can be 
found in Steidel et al. (2003; 2004) and Erb et al. (2005, 
in preparation). 

The QSO spectra we analyzed were obtained between 
1996 and 2003 with either the HIRES echelle spectro- 
graph (Vogt et al. 1994) or the ESI echellette (Sheinis et 
al. 2000) on the Keck I and II telescopes. The HIRES 
spectra were reduced using the standard procedures of 
Tom Barlow's Makee package. The ESI spectra were re- 
duced with the "Dukee" suite of custom IRAF scripts 
and C programs written by KLA and M. Hunt. The 
shape of the continuum was removed from the QSO spec- 
tra interactively by fitting a high-order (~ 30) b-splinc 
to regions of the spectra that appeared to be free from 
absorption. CIV absorption systems in the QSO spec- 
tra were identified interactively by scanning for absorp- 
tion doublets with the right wavelength spacing. Col- 
umn densities were estimated by fitting each doublet's 
absorption profile with two Gaussians in optical depth. 
This is equivalent to Voigt-profile fitting at the low (i.e., 
undamped) column-densities of interest to us. Occasion- 
ally the velocity profiles were unresolved in our spectra. 
If this happened and the doublet ratio suggested that 
the lines were saturated, we obtained a crude estimate 
of column density from the weaker line only. Otherwise 
both lines were used. See Adelberger et al. (2003) for 
further details. Our ability to detect CIV systems de- 
pended on the quality of the QSO spectrum. Table |3 
lists the minimum detected CIV column density in each 
QSO's spectrum. This should be roughly equal to the 
spectrum's CIV detection limit. 

Figure [3] shows the redshift distributions of the differ- 
ent objects in our sample. 

3. GALAXIES' GASEOUS ENVELOPES 

Numerous authors have noticed that the interstellar 
absorption lines in high-redshift galaxies' spectra tend 
to be blueshifted, presumably because they are produced 
by outflowing gas (e.g., Franx et al. 1997; Pettini et al. 
1998, 2001; Frye, Broadhurst, & Benftez 2002; Shapley 
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Table 2. Observed QSOs 



QSO 


Field 


a(2000) 


.5(2000) 


z 


G A B a 


n civ b 


log(AT3? T \n 
° V d V 


PKS0201+113 d 


Q0201 


02 ,l 03 m 46.7 3 


11°34'45" 


3.610 


20.1 





- 


LBQS0256-0000 


Q0256 


02 ,l 59 m 05.6 3 


00°11'22" 


3.324 


18.0 


32 


12.02 


LBQS0302-0019 


Q0302 


03 h 04 m 49.9 3 


-00°08'13" 


3.267 


17.8 


16 


12.03 


FBQS J0933+2845 


Q0933 


09 h 33 m 37.3 3 


28°45'32" 


3.401 


18.8 


34 


12.00 


Q1305kkc9 


Q1305 


13 h 07 m 42.7 3 


29°19'36" 


2.462 


21.3 


7 


13.56 


Q1305kkcl6 


Q1305 


13 h 08 m 06.1 3 


29°22'39" 


2.526 


20.9 


4 


13.03 


Q1305kkc22 


Q1305 


13 h 08 m 11.9 3 


29°25'13" 


2.979 


19.0 


16 


13.02 


Q1422+2309 


Q1422 


14' l 24 m 38.1 s 


22°56'01" 


3.515 


16.5 


74 


11.81 


Q1422+2309b 


Q1422 


14 fc 24 m 40.6 s 


22°55'43" 


3.084 


23.4 


3 


13.45 


Q1623kp76 


Q1623 


16' l 25 m 48.1 s 


26°44'33" 


2.246 


18.5 


13 


12.72 


Q1623kp77 


Q1623 


16 h 25 m 48.8 3 


26°46'59" 


2.529 


16.5 


11 


12.03 


Q1623kp78 


Q1623 


16 h 25 m 57.4 3 


26°44'49" 


2.578 


19.0 


47 


12.14 


Q1623kp79 


Q1623 


16 h 26 m 06.2 3 


26°50'33" 


2.180 


19.2 


12 


12.99 


Q1623BX234 


Q1623 


16' l 25 m 33.7 s 


26°53'44" 


2.464 


19.5 


4 


14.06 


Q1623BX603 


Q1623 


16 h 26 m 02.4 3 


26°43'39" 


2.529 


20.5 


8 


13.32 


HS1700+6416 


Q1700 


17 h 01 m 00.6 3 


64°12'09" 


2.717 


16.0 


57 


11.71 


Q2233+136 


Q2233 


22 h 36 m 27.2 3 


13°57'13" 


3.209 


18.7 


25 


12.09 


Q2342+125 


Q2343 


23 ,l 45 m 22.8 3 


12°45'47" 


2.497 


18.2 


9 


12.57 


Q2343+125 


Q2343 


23 h 46 m 28.3 3 


12°48'58" 


2.579 


17.0 


38 


12.25 


Q2344+125 


Q2343 


23 h 46 m 46.1 3 


12°45'27" 


2.785 


17.5 


11 


12.33 


Q2343BX415 


Q2343 


23 h 46 m 25.4 3 


12°47'44" 


2.578 


20.3 


19 


13.22 


Q2345+000 


Q2346 


23 h 48 m 25.4 3 


00°20'41" 


2.652 


19.7 


3 


13.31 


SSA22dl3 


SSA22a 


22 h 17 m 22.3 s 


00°16'41" 


3.343 


21.6 


6 


13.22 



a QSO magnitude; numbers arc unreliable brighter than our saturation magnitude G ~ 18 
b Number of detected CIV systems 

c Log of the minimum column density (in cm -2 ) among the detected CIV systems 
d CIV absorber catalog not constructed owing to gaps in our echelle spectrum 
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Fig. 3. — Redshift distributions for different objects in our 
sample. 



et al. 2003; Erb et al. 2004). Little is known, however, 
about how far the outflowing gas might propagate. This 
section shows that the absorbing material often lies at 
radii approaching 80 proper kpc and that weaker lines 
can be detected much farther away. It does not clearly 
show, however, that the absorbing material detected at 
large radii is part of an outflow. Other scenarios are 
considered in §§0]and|Hl 



3.1. Radii r ^ 80 proper kpc 

A ~ 40 kpc radius for galaxies' strongly absorbing 
material is implied by the presence of foreground galax- 
ies' absorption lines in the spectra of nearby background 
galaxies. The left panel of Figure 0] illustrates the idea. 
The right panel shows one example from our survey. Two 
galaxies with redshifts z = 1.60, z = 2.17 and angular 
separation 8 = 2" were placed on the same slit on one of 
our multislit masks. Light from the background galaxy 
passed within 17 proper kpc of the foreground galaxy as 
it traveled towards earth, and this was close enough for 
the foreground galaxy's absorption lines to be imprinted 
in its spectrum. The absorbing material must therefore 
have a radius of at least 17 kpc. 

Similar absorption lines are observed in galaxy pairs 
with impact parameters of up to 40 kpc, but the absorp- 
tion weakens significantly at greater separations. This 
can be seen by measuring the dependence on angular 
separation 9 of the the mean absorption produced by 
foreground galaxies. We calculated the strength of the 
mean absorption in two bins of angular separation. For 
each bin, we identified all the galaxy pairs in our sample 
(18 for 1" < 9 < 5", 64 for 5" < 8 < 10"; 10" cor- 
responds to 81 proper kpc at z = 2.5 for a cosmology 
with — 0.3, Qa = 0.7, h = 0.7), collected the back- 
ground spectra, masked parts of the spectra that were 
contaminated by sky lines or by any of the background 
galaxy's ten strongest spectral features 4 , shifted each 
spectrum into the rest-frame of its foreground galaxy, 

4 Ly/3/OVI, Lya, Sill A1260, OI A1302, CII A1335, SilV A1398, 
Sill A1527, CIV A1549, Fell A1608, A1II A1671 
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Fig. 4. — Left: Schematic view of galaxy/galaxy absorption. 
A foreground galaxy, B, is surrounded by a gaseous envelope that 
may be expanding. Light from the background galaxy, A, passes 
through this envelope and is imprinted with absorption lines. (A 
may have an expanding envelope of its own but this is suppressed 
for simplicity.) Right: Example LRIS spectra of a galaxy pair with 
2" separation and redshifts z = 1.61, 2.17. Spatial position varies 
along the x axis and wavelength increases towards higher y. The 
Lyman-a emission line of the higher redshift object is visible at 
the bottom of the rightmost spectrum. The gas associated with 
the lower redshift object, on the left, produces absorption lines 
(labeled) in the spectra of both objects. Its radius must therefore 
exceed 17 kpc (i.e., 2" at z = 1.61). 
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Fig. 5. — Mean gas absorption as a function of galactic im- 
pact parameter. Top panel: mean spectrum of the foreground 
galaxies in pairs with angular separation 1" < 8 < 5" (8 ^ b ;$ 40 
proper kpc). Their absorption lines are sensitive to the gas den- 
sity at impact parameters smaller than their half-light radii, i.e., 
8 ^ 0.2". The names of different absorption and emission lines are 
indicated. The darker shaded region at the bottom of the panel 
is the bootstrap error spectrum, estimated by recalculating the 
mean spectrum for random ensembles of foreground galaxy spec- 
tra. Middle panel: mean absorption observed in the spectrum of 
the background galaxy at the redshift of the foreground galaxy. 
This panel is noisier, owing to our masking of interstellar features 
and (at bluer wavelengths) to Lyman-a forest absorption in the 
spectra of the background galaxies. Apparent features should be 
viewed skeptically. Comparison to the error spectrum shows that 
some of them are not significant (e.g., absorption near 1375A, emis- 
sion near 1580 A). Bottom panel: similar to the middle panel, but 
this is the mean absorption produced in the spectra of background 
galaxies with angular separation 5" < 8 < 10" (40 ,$ b ,$ 80 proper 
kpc). 



spline-interpolated the rest-frame spectra onto common 
abscissae, removed the shape of the continuum by forcing 
each spectrum to have a constant 140A running-average 
flux, and averaged the spectra together with inverse- 
variance weighting after first rejecting the highest and 
lowest 10% of the flux values at each output wavelength. 
Finally we restored a realistic continuum shape by mul- 
tiplying the resulting spectrum by the mean running- 
average continuum shape of the sample. 

Figure [S] shows the result. The strength of the inter- 
stellar CIV absorption line appears roughly constant for 
9 ^ 5" (~ 40 proper kpc at z = 2.2) but drops sharply af- 
terwards. 5 Chen, Lanzetta, & Webb (2001) found similar 
results at lower redshifts. It would be wrong to conclude 
from Figure however, that the observed CIV absorp- 
tion at 5" < 8 < 10" is weak to the point of insignifi- 
cance. The estimated rest-frame equivalent width of the 
absorption, W (doublet) ~ 0.7A (i.e., Wyi548) > 0.4A, 
or Aciv ~ 10 14 cm -2 ), is large by intergalactic stan- 
dards. According to Steidel (1990), only 1.2 ± 0.4 CIV 
systems with W (1548) > 0.4A are found per unit red- 
shift at 2 z <, 3. Since our color-selection criteria find 

5 Note that part of the CIV feature in the top panel of Figure l5l 
is stellar P-Cygni absorption. 



9 galaxies per arcmin 2 per unit z at the same redshifts 
(e.g., Steidel et al. 2004), they would account for the ma- 
jority of all CIV-systems with W (1548) > 0.4A if each 
galaxy produced absorption with Wo (1548) > 0.4A out 
to 9 = 10" = 0.167': 9 x tt(0.167') 2 ^ 0.8 CIV absorbers 
per unit redshift. In fact, as we will see below (Figure^, 
only about half of the galaxies within 10" (~ 80 kpc) are 
associated with CIV absorption so strong. This implies 
that the galaxies in our sample can account for roughly 
one-third of CIV lines with W (1548) > 0.4A. Galaxies 
missed by our UV-color-selection techniques could easily 
account for the remainder. 

The large equivalent width Wo (doublet) ~ 2.7A of 
the absorption at 9 < 5" has an interesting implica- 
tion. CIV can attain so high an equivalent width only 
if peculiar velocities spread the absorption over a large 
range of wavelengths. The minimum velocity spread is 
At) = 260 km s , but this assumes a contrived situa- 
tion where each line in the doublet is saturated and has 
a boxcar (ie, maximally compact) absorption profile. In 
realistic situations the CIV absorption will have compli- 
cated substructure (e.g., Pettini ct al. 2002) and a sig- 
nificantly larger velocity spread will be required to pro- 
duce Wo (doublet) ~ 2.7A. Our lower limit to the velocity 
spread (Av — 260 km s _1 ) is smaller than the likely es- 
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cape velocity at 40 kpc from a galaxy with M ~ 1O 12 M 
(e.g., Adelberger et al. 2004), but it nevertheless shows 
that significant chaotic motions in the gas extend to large 
radii. 

Further evidence comes from one of the two cases in our 
survey where a foreground galaxy lies within 15" of the 
QSO sightline and the QSO is lensed into two detectable 
components with ~ 1" separation. In this case, shown 
in the left panel of Figure the galaxy's gas produces 
significantly different absorption profiles in the spectra 
of the two QSO components. The implied substructure 
on half-kpc scales would be erased by thermal motions 
in a few dozen Myr if it were not maintained somehow. 
The gas along the two sightlines might have supersonic 
relative motions or might be gravitationally confined in 
two separate mini-halos. Substructure of this sort is rare 
in the low-density intergalactic medium (Rauch et al. 
2001b) but common among CIV systems (Rauch, Sar- 
gent, & Barlow 2001a). 

3.2. Radii r ^ 80 proper kpc 

Although the metal absorption lines continue to 
weaken beyond r ~ 80 kpc, they are not always ab- 
sent. High signal-to-noise QSO spectra reveal that ab- 
sorbing gas extends out to 100-200 proper kpc in at least 
some cases. Figure [7| shows three. Absorption from 
metal-enriched, multiphase gas is detected in each. (See 
§§ 2.5.13 and 2.5.15 of Simcoe et al. 2002 for a discus- 
sion of the physical conditions in two of these absorption 
systems.) 

These three cases are not typical, but nor are they ex- 
traordinarily rare. Figure|S]shows the total detected CIV 
column density within 100-500 km s _1 of the 29 galaxies 
in our NIRSPEC sample that lie within 1.0/i -1 projected 
comoving Mpc of a QSO with detected CIV 6 . The CIV 
column densities of the three examples in Figure {7\ arc 
the second, third, and seventh highest in this 29 object 
sample. 

One way to show that the detected metals are asso- 
ciated with a nearby galaxy, and are not chance align- 
ments, is to measure the galaxy-CIV correlation strength 
£ gc on small spatial scales. As pointed out by Adelberger 

et al. (2003), £ gc cannot easily exceed fgg 2 fc/ 2 , the ge- 
ometric mean of the galaxy-galaxy and CIV-CIV clus- 
tering strengths, unless the existence of galaxies affects 
(or is affected by) the presence of nearby CIV systems. 
Since £ gg and £ cc are roughly equal (Quashnock & van- 
den Berk 1998; Adelberger et al. 2003), a value of £ gc /£gg 
significantly greater than unity on small scales would im- 
ply a direct relationship between the galaxies and CIV 
systems. 

We estimated £ gc /£g g with the following approach. For 
each CIV system, we counted the observed number of 
NIRSPEC galaxies with velocity separation \Av\ < v max 
and transverse separation i? m i n < R < -Rmax- Call 
these galaxies the CIV system's neighbors. (We took 
Vmax = 500 km s _1 and i? max ~ Ih^ 1 comoving Mpc, 
since these are roughly the maximum separations ex- 
pected for a galaxy's ejecta.) Given the galaxies' known 
angular distances to the QSO, the expected number of 

6 A small fraction of our QSOs were too faint to allow us to 
detect significant numbers of CIV systems. 



neighbors is (see, e.g., Adelberger 2005) 
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where the sum runs over all NIRSPEC galaxies in the 
QSO's field, Az is the redshift difference corresponding 
to v max , P( z ) is the survey selection function normal- 
ized so that J °° = 1, £ is the correlation func- 
tion, and Yij is the distance between the CIV-system at 
redshift Zj and a point at redshift z with the galaxy's 
angular separation 8j. The number of neighbors ex- 
pected in the absence of clustering, n 6 , can be calcu- 
lated from equation 0] with £ = 0, allowing us to_esti- 
mate the galaxy-CIV clustering strength through £ gc = 
— 1 + J2i n i/ J2i ™i- This is the mean value of £ gc within 
a complicated volume consisting of many parallel skew- 
ers of length 2w max . To estimate the mean value of the 
galaxy-galaxy correlation function within the same vol- 
ume, we inserted £ = £ gg = (r/4/i comoving Mpc) -1 ' 6 
(e.g., Adelberger et al. 2005) into equation^ integrated 
numerically, called the result rip, and used the formula 

£gg = — 1 + ^f S / Si ■ The estimates of £ g c and £ gg 
do not depend on any assumptions about the redshift 
selection function of the CIV systems or the angular se- 
lection function of the galaxies. This is fortunate since 
neither is well known. 

Figure shows £ gc /£g g for ^min,-Rmax = 0,0.4 and 
0.4, 1.0/i -1 comoving Mpc. Our data are consistent with 
no increase in £ gc over £ gg on small scales. They are 
therefore consistent with the idea that galaxies and CIV- 
systems are spatially correlated only because they trace 
the same large-scale structure. We cannot rule out a 
more direct association, however, both because the error 
bars in Figure are large and because the statistical test 
itself can detect only the strongest direct associations of 
galaxies and CI V-sy stems. 

We conclude this subsection by emphasizing one of 
its corollaries: many of the strongest intergalactic ab- 
sorption lines are produced by material that lies close 
to star- forming galaxies. Figure ITUI helps illustrate the 
point. The figure's vertical lines mark the redshifts 
of strong CIV and OVI absorption in the spectrum of 
HS1700+6416. 7 Circles mark the positions of galaxies 
near the QSO sightline. Small regions with a wealth of 
galaxies and intergalactic metals tend to be surrounded 
by vacant expanses with neither. The correspondence 
between star-forming galaxies and strong metal lines is 
not perfect, but one could do worse than assume that the 
presence of one implies the presence of the other. 

4. THE GALAXY-CIV CROSS-CORRELATION FUNCTION 

The previous section discussed the relationship be- 
tween galaxies and metal absorption lines for a small 
number of specific cases. This section takes a more sys- 
tematic approach, calculating the cross-correlation func- 
tion between the positions of galaxies and CIV absorp- 
tion systems. 

4.1. Large scales 



7 The redshifts of OVI absorption were taken from Simcoe et al. 
(2002); none of the other QSOs in our sample were analyzed by 
them, and we have not attempted to find OVI systems ourselves. 
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Fig. 6. — Substructure in metal-enriched gas near galaxies. Each panel is for a single galaxy that lies within 15" of a background QSO 
that is lensed into two components with ~ 1" separation. The proper projected distance from galaxy to QSO sightline (b) and proper 
projected separation of the two QSO components (Ar) are indicated above the panels; Ar is taken from Rauch et al. 2001a. The ±1(T 
range of the galaxy's redshift is indicated by the rectangle at the bottom of each figure. The Lyman-a and metal-line absorption near the 
galaxy's redshift is shown in the stacked spectral segments, which are drawn shaded for one QSO component and dashed for the other. 
Significant differences are apparent in the left-hand panel but not the right. (In the right hand panel, the Lyman-a forest falls in the same 
part of the QSO spectrum as the SilV and CII absorption lines, making it difficult to determine if these features are present or absent.) 
The recently obtained galaxy redshifts near Q0142-100 have been used only for this plot, and this field has consequently been omitted from 
Tables and El 



Our approach on large spatial scales is similar to that of 
Adelberger (2005): we assume that the cross-correlation 
function has the form £(r) = (r/ro) -1 ' 6 (e.g., Adelberger 
et al. 2003); count the number n o bs(0,£) of galaxy-CIV 
pairs in our sample that have A9 < 300", \AZ\ < I for 
I — 20, 40/i -1 comoving Mpc 8 ; calculate the expectation 
value of n o bs(0, 20)/n o bs(0, 40) as a function of Tq with 

8 AZ is the comoving separation in the redshift direction be- 
tween the galaxy and CIV system; i.e., it is the comoving distance 
between the galaxy's redshift z g and the CIV system's z c . 



the equation (Adelberger et al. 2005; Adelberger 2005) 
n obs (0,1) \ _ EuJpdZll+m^Z)} 

^• M )/~E, J J? < ^i + {(%2)]' 

where £(R,Z) = i[(R 2 + Z 2 ) 1 /\R l] = {l+z i )D A {z i )e ij , 
Da(z) is the angular diameter distance, and the summa- 
tion is over all galaxies (i) and CIV-systems (j); and fi- 
nally take as our best estimate of ro the value that makes 
the right-hand side of equation [S] equal the observed ra- 
tio n o bs(0, 20)/n o bs(0, 40). As discussed in Adelberger et 
al. (2005) and Adelberger (2005), the resulting estimate 
of ro is insensitive to angular selection effects and only 
weakly dependent on the assumed shape of the selection 
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Fig. 7. — Absorption lines in QSO spectra at the redshifts of nearby galaxies. Each of the three panels (left, middle, right) is for a single 
galaxy, with proper impact parameter indicated above. The galaxy's redshift, measured with NIRSPEC, is marked by the short bar in the 
center of the panel. The width of the bar, ±60 km s , is equal to the redshift's la uncertainty. The redshifts of neighboring galaxies 
are indicated with the other bars. The neighboring galaxies were not observed with NIRSPEC and so their redshifts are more uncertain, 
with ±lcr limits roughly 5 times broader than the bars. The bars' heights are proportional to the galaxies' impact parameters. The HI 
absorption near the galaxy's redshift is shown at the bottom of the panel. Accompanying metal-line absorption, if any, is shown above. 
Redshifts of detected metal absorption are marked with downward pointing arrows near the top of the panel. The velocity range shown in 
each panel, roughly ±900 km s — 1 , is small enough that absorption lines visible anywhere in a panel could conceivably be produced by gas 
flowing out of the galaxy in the middle. Some of the lines (e.g., OVI, CII) fall within the Lyman-a forest, and observed absorption can be 
attributed to them with confidence only if other metal-lines are present at the same redshift. 



function — two significant advantages given the irregular 
selection criteria we employed (§[2}. 

Figure ITTI shows the estimated cross-correlation length 
as a function of the CIV system's column density. To 
estimate the uncertainty, we broke our galaxy and CIV 
catalogs into many smaller sub-catalogs by rejecting a 
random fraction p = 0.5-0.8 of the sources, observed 
how the dispersion in best-fit rn among the sub-catalogs 
depended on the number of objects in the catalogs, and 
extrapolated to the full catalog sizes. This accounts for 
random uncertainties but not completely for cosmic vari- 
ance. 



Also drawn in Figure ^2 is the la confidence inter- 
val for the galaxy-galaxy correlation length at similar 
redshifts, taken from Adelberger et al. (2005). The 
good agreement of the galaxy-galaxy and galaxy-CIV 
correlation lengths for Nciv ~ 10 125 cm" 2 shows that 
bright star-forming galaxies and stronger CIV systems 
have similar spatial distributions, reinforcing the idea 
(§ that they are often the same objects. Adelberger 
et al. (2003) reached a similar conclusion. 

4.2. Small Scales 

As mentioned above, an observed spatial association 
of galaxies and CIV systems does not by itself imply 



HIGH-REDSHIFT GALAXIES AND INTERGALACTIC ABSORPTION LINES 



9 



10" 



0.2 0.4 0.6 0.8 
b / comoving h" 1 Mpc 



20 



|Av|<500 | | 
|Av|<200 
|Av|< 100 [J " 



Jo 



10 11 10 12 10 13 10 14 10 15 10' 6 
N'°{ / cm" 2 



Fig. 8. — Left panel: Column density of the strongest detected 
CIV system within ±200 km s — 1 of the galaxy as a function of 
galaxy impact parameter b. Large points show detections. Small 
points with arrows show upper limits. We assumed that the upper 
limit would be equal to the column density of the weakest detected 
CIV system in the QSO's spectrum. Right panel: Distribution 
of the total detected CIV column density near galaxies that lie 
within l/i -1 comoving Mpc of the QSO sightline. The total column 
density was calculated by summing the column densities of all CIV 
systems within 100, 200, or 500 km s — 1 of the galaxy's redshift. 
For simplicity, the non-detections are all placed in the bin with 
JVciv < 10 12 cm -2 . This panel provides some indication of how 
the left panel would change for velocity limits other than ±200 km 
s — L The same galaxies were used in both panels; all have precise 
(i.e., near-IR nebular) redshifts. 
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Fig. 10. — Redshifts of galaxies near the sightline to QSO 
HS1700±6416 compared to the redshifts of CIV and OVI absorp- 
tion systems in the QSO's spectrum. Circles mark galaxy redshifts 
and impact parameters. The area of each circle is proportional 
to the galaxy's apparent luminosity in the 1Z band. Vertical lines 
mark the redshifts of OVI and CIV absorption. Thicker lines in 
the CIV panel mark systems with multiple components. OVI ab- 
sorption redshifts are taken from Simcoe, Sargent, & Rauch (2002). 
The absorption line catalogs are reasonably complete down to their 
limiting equivalent widths, but, owing to imperfect color-selection 
criteria and limited time for follow-up spectroscopy, it is unlikely 
that the galaxy catalogs contain more than half of the galaxies 
brighter than 1Z = 25.5 at these redshifts. (We obtained redshifts 
for 15 of the 20 photometric candidates within 60" of the sightline 
to this QSO, so most of the expected incompleteness is due to our 
color-selection criteria. See, e.g., Adelberger et al. 2004.) 
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Fig. 9. — Estimated ratio of mean galaxy-CIV and galaxy-galaxy 
clustering strength as a function of galaxy impact parameter b. 
The correlation functions are averaged over ±500 km s — 1 for each 
galaxy. Points show the ratios for CIV systems in two column- 
density ranges, as indicated. The lcr error bars were calculated 
under the assumption that the variance in the number of pairs N 
is Var(AQ ~ 1.56JV 1 ' 24 , an approximation that Adelberger (2005) 
found to be accurate for galaxy pair counts in the Lyman-break 
survey. 
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Fig. 11. — Galaxy-CIV cross-correlation length ro as a function 
of CIV column density. Circles show the measurements of sec- 
tion ^] The number near each point indicates the number of CIV 
systems with column densities in the bin's range. The shaded re- 
gion in the background shows the galaxy-galaxy correlation length 
ro appropriate to this sample (Adelberger et al. 2005). 
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that the detected metals are part of a galactic super- 
wind. Many other processes could place metals near 
galaxies. One way to learn about the relationship be- 
tween intergalactic CIV and nearby galaxies is to look 
for anisotropies on small scales in the redshift-space 
galaxy-CIV correlation function. Since these anisotropies 
are produced by peculiar velocities, they should reveal 
whether the CIV tends to be falling towards, flowing 
away from, or orbiting around the nearest galaxy. Fig- 
ure El shows the observed spatial separation of every 
galaxy-CIV pair in our NIRSPEC sample that had co- 
moving redshift separation \AZ\ < 10ft _1 Mpc and co- 
moving angular separation b < lh^ 1 Mpc. Vertical lines 
mark the impact parameters of each galaxy-QSO pair. 
Circles on these lines show the velocities (relative to the 
nearby galaxy) of any detected CIV absorption in the 
QSO's spectrum. Only galaxies with NIRSPEC redshifts 
are shown in the plot, since precise redshifts are crucial 
to our arguments. 

The peculiar velocities of CIV systems relative to 
galaxies can be estimated roughly from this plot as fol- 
lows. The expected number of galaxy-CIV pairs with 
impact parameter b and radial separation \AZ\ < Z\ is 
proportional to [1 + £(&, < Zx)]Zi, where 



i(b,<Z 1 ) 



1 



dZ£(&, Z) 



(6) 



and £(&, Z) is the galaxy-CIV correlation function. On 
average, then, half of the plot's galaxies with impact 
parameter b should have comoving redshift separation 
\AZ\ < Z 1/2 (b) with 



Z 1/2 (b) 



l + < z max ) 
2 i + i(b,<z 1/2 ) 



(7) 



where \AZ\ < Z max is the range of redshift separations 
shown on the plot. If peculiar velocities were negligi- 
ble, £(6, Z) would be equal to the real-space galaxy-CIV 
correlation function £ r , which must be isotropic in an 
isotropic universe: £{b,Z) — C[{b 2 + Z 2 ) 1 / 2 ]. Since 
£ r (r) = (3.7/i _1 Mpc/r) 16 is a reasonably good approx- 
imation to the real-space correlation function (see Fig- 
urc|51), one can solve for Zi/ 2 {b) numerically by inserting 
this correlation function into equations |7| and 

The dotted triangular envelope in Figure El shows 
Zi/ 2 (b). Half of the CIV systems would be expected to 
lie within the dotted envelope in the absence of peculiar 
velocities. In fact 66 of 81 lie outside. Peculiar velocities 
appear to be substantial. The observations seem incon- 
sistent with an infall model, since infall tends to compress 
correlation functions in the redshift direction; it would 
place more than half of the CIV systems inside the en- 
velope. The random motions of orbiting clouds would 
move some CIV systems outside the envelope through 
the finger-of-god effect, but it is unclear if the expected 
size of the displacements (±1.4/i _1 comoving Mpc for 
aid = 140 km s _1 ; e.g., Adelberger et al. 2005) is large 
enough to account for our observations. The small sam- 
ple size prevents us from drawing firm conclusions, but 
outflows moving at several hundred km s _1 might be- 
come the favored explanation if the correlation function 
remains so anisotropic as the sample size grows. 

5. INTERGALACTIC HI NEAR GALAXIES 



After noticing that none of the three galaxies closest to 
background QSOs in their sample seemed to be associ- 
ated with strong Lyman-a absorption lines, Adelberger 
et al. (2003) speculated that superwinds extending to 
~ 0.5/i _1 Mpc might drive most intergalactic hydrogen 
away from young galaxies. Our larger sample puts this 
speculation to rest. The top panel of FigurelTSlshows the 
measured transmissivity of every Lyman-a forest pixel 9 
in our QSO spectra that lies within l/i -1 comoving Mpc 
of a galaxy with precisely known (i.e., near-IR nebular) 
redshift. The distances shown on the x-axis are impre- 
cise, because they were calculated from redshift differ- 
ences under the assumption that peculiar velocities are 
negligible, but nevertheless it is clear that strong HI ab- 
sorption is common near the galaxies. 

The bottom panel of Figure El shows the same data 
with the pixels grouped by galaxy. Each point marks 
/iMpci the mean Lyman-a transmissivity of the pixels 
within l/i -1 comoving Mpc of a single galaxy, and the 
error bar shows the r.m.s. spread among these pixel 
values. Table |3 lists the impact parameters and red- 
shifts for the galaxies in this figure. Also listed is c, 
the smoothed mean transmissivity in the QSO's spec- 
trum at the galaxy's redshift, calculated by convolving 
the continuum-normalized QSO spectrum by a boxcar 
with width 100A and fitting a 2nd-order polynomial to 
the result. 

Figure ITU shows that the mean transmissivity declines 
monotonically as one approaches a galaxy. This appears 
inconsistent with the earlier result of Adelberger et al. 
(2003; dashed circles). The present analysis differs from 
theirs in three main ways: our redshifts are more precise, 
our sample is larger, and our galaxies' typical redshift is 
lower. Each difference could contribute to the change 
in this plot. The good agreement of the results for the 
NIRSPEC-only and full samples in Figure El suggests 
that the redshifts' precision might not be primarily re- 
sponsible for the change. This suspicion is reinforced by 
our NIRSPEC observations of galaxies in the sample of 
Adelberger et al. (2003). Although two of their three 
galaxies with r < 0.5ft -1 comoving Mpc had redshifts 
that placed the nebular lines outside of atmospheric win- 
dows, the near-IR nebular redshift we obtained for the 
third confirmed the reported lack of nearby Lyman-a 
absorption. 10 Another possibility is that genuine evolu- 
tion from z ~ 3 to z ~ 2 is partly responsible for the 
increased absorption near galaxies in this sample. The 
thought is not absurd: at lower redshifts, galaxies have 
marginally slower outflows (see Figure ^) and a larger 
fraction of intergalactic HI absorption is produced by 
gas in deep potential wells that is more difficult to dis- 
rupt. These changes are not enormous, however, and it 
seems unlikely to us that they would have a dominant 
effect. The most obvious difference from Adelberger et 
al. (2003) is the sample size, especially at small sepa- 
rations. The following simple argument suggests that a 
statistical fluctuation could plausibly be responsible for 
lack of absorption at r < O-5/i" 1 comoving Mpc reported 

9 We excluded pixels within damped Lyman-a systems, within 
Az = 0.05 of the QSO, or within parts of the Lyman-a forest that 
could be contaminated by Lyman-/3 absorption from gas at higher 
redshifts. 

10 though it altered the galaxy's redshift by more than antici- 
pated by Adelberger et al. (2003). 
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Fig. 12. — Spatial separations of galaxy-CIV pairs in the angular (ft) and redshift (AZ) direction. Vertical lines show the projected 
distance between each galaxy in the NIRSPEC sample and its nearest QSO. The QSO spectra allow us to detect any CIV-absorbing gas 
along these lines. Circles show the actual locations of detected CIV absorption, with one circle at the redshift of each resolved absorption 
component. In the absence of peculiar velocities, the distribution of CIV systems around galaxies would be spherically symmetric on 
average and (as a result) half of the circles would be expected to lie inside the dotted triangular envelope (see §BJ. In fact 65 of 81 lie 
outside. The uncertainty in each galaxy-CIV pair's redshift separation (i.e., in each circle's vertical position) is too small to account for 
this (±1<7 size indicated by two-headed arrow labeled Az). Peculiar velocities must therefore be substantial. If the CIV systems were 
orbiting in the galaxy's halo, their observed redshift separations would be displaced relative to their true separations by a random amount 
comparable to their galaxy's velocity dispersion <r ltJ . The approximate size of ±<Tid (140 km s^ 1 ) is indicated by the labeled double-headed 
arrow. If the CIV were infalling, the CIV systems' observed (i.e., redshift space) positions would be artificially compressed towards the x 
axis by an amount comparable to aid- If it were flowing outwards, the observed positions would be displaced away from the x axis by an 
amount that depends on the wind speed. The two-headed arrow labeled ^ w i n d shows the approximate size of ±500 km s — 1 . 



by Adelberger et al. (2003). Suppose that half of the 
galaxies with r < 0.5/i _1 comoving Mpc are heavily ob- 
scured (/ = 0) and half are unobscured (/ = 1). This is a 
rough approximation to the data in the bottom panel of 
Figure lTSl Then the mean and standard deviation of the 
average transmissivity for a three-object sample would 
be fs = 0.5 ± 0.29. Although an observation of / 3 = 1 
superficially seems inconsistent with this_scenario, differ- 
ing by almost 2a from the expectation — 0.5, in fact 
it should occur one time out of eight. The real world is 
only vaguely similar to our simple model, but the point 



is that large fluctuations occur often in samples drawn 
from a bimodal underlying distribution. The results in 
this paper and in Adelberger et al. (2003) are clearly not 
the same, but they are less inconsistent than one might 
naively suppose. 

In any case, although strong HI absorption is the norm 
at small separations, we continue to find some galaxies 
with little absorption. The existence of these galaxies 
is surprising: one would expect all galaxies to reside in 
dense regions with significant hydrogen and an elevated 
neutral fraction from the p 2 dependence of the recom- 
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Table 3. Galaxies with impact parameter b < comoving Mpc a 
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a Unobserved quantities are marked "-" ; quantities we could not measure owing to noise in the data are marked "x." 

b Comoving impact parameter, h~ 1 Mpc 

c Redshift 

d Mean transmissivity of pixels within lh~ 1 comoving Mpc 

e Mean transmissivity of random pixels with similar redshifts in the QSO's spectrum 
f Total detected CIV column density within 200 km s~ 1 of the galaxy redshift, cm -2 

g Velocity difference between interstellar absorption lines and nebular emission lines, km s . Typical uncertainty ~ 200 km s _1 . 
h Velocity difference between Lyman-a and nebular emission lines, km s — 1 . Typical uncertainty ~ 100 km s _1 . 
'Stellar mass, 10 10 AfQ. Typical uncertainty ~ 0.15 dex 
JStar formation rate, Mq yr -1 . Typical uncertainty ~ 0.2 dex 
k Dust reddening E(B — V). Typical uncertainty ~ 0.03 
'Age, 10 s yr. Typical uncertainty ~ 0.3 dex 

"Time constant of exponentially declining star-formation history, 10 s yr. Typical uncertainty ~ 1 dex 



bination rate (e.g., Croft et al. 2002; Kollmeier et al. 
2003; Maselli et al. 2004; Desjacques et al. 2004). 
Figure El illustrates this point by comparing the ob- 
served distribution of HI absorption among the galax- 
ies in Figure 1131 to a preliminary prediction from the 
windless smooth-particle hydrodynamical (SPH) simu- 
lation of Kollmeier et al. (2005, private communica- 
tion; this simulation is similar to that of Kollmeier et 
al. 2003 but is matched to the mean redshift of this 
sample). If the real world resembled the simulations, a 
24 galaxy sample would have n Q b s = 7 or more mem- 
bers with decrement D < 0.2 only two times in a 
thousand (i.e., X^ ba=7 exp(-n CX p)n™°£ a /n obs ! ~ 0.002, 
where n eX p — 1-75 is the predicted mean number of galax- 
ies with D < 0.2 in a 24 object sample). 

Although a number of explanations are possible, the 
data in Figure HTl seem qualitatively consistent with the 
idea that the galaxies are surrounded by anisotropically 
expanding winds. Leaking out along the paths of least 
resistance, these winds would turn moderate absorption 



into weak absorption while leaving the strongest absorp- 
tion systems in place. 

Ionizing radiation from any neighboring QSOs might 
also preferentially eliminate weaker absorption systems. 
As Figure 1161 shows, many of the galaxies with /iMpc 
have bright QSOs in their vicinity. Although these QSOs 
have similar or larger redshifts than the galaxies, as is 
required for their (presumably) shortlived radiation to 
have an observable effect on the galaxies' surroundings 
(e.g., Figure 12 of Adelberger 2004), the prevalence of 
neighboring QSOs is not significantly different for galax- 
ies with /iMpc > 0.5 and /iMpc < 0.5. Four galaxies of 9 
with /iMpc > 0.5 have a QSO within Av = 3000 km s _1 , 
compared to 2 of 15 with /im P c < 0.5. The one-sided P- 
valuc from Fisher's (1934) exact test is 0.11, and would 
increase (i.e., become less significant) if we changed Av 
to 2000 km s" 1 (P = 0.25) or 4000 km s" 1 (P = 0.21). 
These numbers could be misleading, however, since even 
a profound QSO influence might be difficult to detect 
with significance in a sample so small. More convincing 
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Fig. 13. — Top panel: observed Lyman-a absorption for every 
QSO-spectrum pixel whose apparent (i.e., redshift space) distance 
to the nearest galaxy is r < lh" 1 eomoving Mpc. The circles 
are larger for lower-resolution spectra; their area is proportional 
to the wavelength size of the pixel and the smallest pixels have 
AA ~ 0.025A. Bottom panel: the same data grouped by galaxy. 
The mean transmissivity for the pixels associated with each galaxy 
is shown on the y axis. Circles are placed on the x axis accord- 
ing to the galaxy's impact parameter. The area of each circle is 
proportional to the galaxy's apparent luminosity in the 1Z band. 
Error bars show the rms spread in flux among the pixels with 
r < l/i -1 Mpc. Darker circles are galaxies from Adelberger ct al. 
(2003), re-analyzed after obtaining more-precise galaxy redshifts 
from near-IR nebular lines. We were able to obtain confirming 
near-IR nebular redshifts for only one of the galaxies from that 
paper with r < 0.5/i -1 eomoving Mpc; the other two had redshifts 
that made near-IR spectroscopy impractical. 
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2 4 6 8 
r / h _1 eomoving Mpc 

Fig. 14. — Top panel: Mean transmissivity as a function of galac- 
tocentric distance. Points at < r < 0.5h~ 1 eomoving Mpc show 
the mean transmissivity of all spectral pixels whose redshift-space 
position lay within 0.5k" 1 eomoving Mpc of a galaxy. To scale to 
a common global transmissivity (/) = 0.765, the transmissivities 
of each pixel were multiplied by (f)/c(z) before averaging, with 
c(z) the redshift-dependent smoothed mean transmissivity in each 
QSO's spectrum (§|oJ. The error bar shows the lc uncertainty, es- 
timated by dividing the observed galaxy-to-galaxy scatter in trans- 
1/2 

missivity by N ^ . The points and error bars at larger separations 
are defined similarly. Darker circles show the result if we restrict 
ourselves to galaxies with near-IR nebular redshifts. Lighter cir- 
cles are for the full sample. For comparison, dotted circles show 
the result of Adelberger et al. (2003) scaled from (/) = 0.67 to 
(/) = 0.765. Bottom panel: the number of galaxies at each sepa- 
ration in each sample. 



is the fact that the QSOs are too faint to alter the ioniza- 
tion balance of material that lies ~ 30 proper Mpc away 
(i.e., At; = 2000 km s^ 1 ; see, e.g., Adelberger 2004). 
Their ionizing radiation could explain the galaxies' weak 
HI absorption only if the QSOs' redshifts differ signifi- 
cantly from our estimates. 

A search for other explanations would have to explore 
possible faults with the data themselves. Precise mea- 
surements of a galaxy's redshift are required to recognize 
that it is associated with a narrow absorption feature in 
a QSOs spectrum. Any errors in the redshift will tend 
to move the galaxy away from the feature, artificially 
reducing the apparent absorption. If such errors were 
responsible for galaxies with weak absorption, one would 
expect those galaxies always to have strong absorption 
features nearby. Figure IT?)! shows the nearby absorption 
features for all nine of the galaxies from Figures[T3land[T5l 
that had transmissivity /iMpc > 0.5. Strong nearby HI 
absorption is present for about half of the galaxies, but 
systematic wavelength errors could make our data consis- 
tent with the windless simulation only if their magnitude 
is much worse than the value a v = 60 km s _1 that we 
estimate from repeated observations of individual galax- 
ies. 

6. CHARACTERISTICS OF POSSIBLE SUPERWIND 
GALAXIES 

The previous sections discussed two characteristics of 
intergalactic absorption that might reflect the influence 
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SPH simulation — no winds 

(Kollmeier et al.) 



Observed 



0.2 0.4 0.6 0.8 1 
flux decrement 

Fig. 15. — Comparison of observed and expected Ly-ct ab- 
sorption near star-forming galaxies. The abscissa shows the flux 
decrement 1 — /iMpci where /iMpci the mean transmissivity among 
all pixels within lh^ 1 eomoving Mpc of th e galaxy, is the quantity 
shown in the bottom panel of Figure [TBI In windless SPH simu- 
lations, almost every sightline near a galaxy has strong Lyman-a 
absorption. This is not observed. 
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Fig. 16. — HI absorption near the redshifts of the 9 sightline galaxies with the weakest absorption in Figure [T^l The lighter shaded 
region shows the HI absorption in the QSO's spectrum. The darker shaded box in the background is centered on the galaxy's redshift. 
Its width marks the range of pixels that lie within lh -1 comoving Mpc of each galaxy. The number to the upper left of each panel gives 
each galaxy's impact parameter in h~ 1 comoving Mpc. The redshifts of neighboring QSOs with Gab < 20.5 in each field are marked with 
stars; the QSOs whose spectra are displayed lie at redshifts too high to fall on this plot. The neighboring QSO in the upper right panel 
and two lowest left panels has Gab ~ 20.5 and lies ~ 9' (4 proper Mpc) from the sightline; the two neighboring QSOs in the upper left 
panel have Gab = 18.5, 19.5 and lie ~ 100", 600" from the sightline. The black square in each panel has a width of 120 km s , twice 
the estimated lcr uncertainty in each galaxy's redshift. The largest positive and negative changes in redshift estimate we have encountered 
when re-observing a galaxy with NIRSPEC, 124 km s _1 and —95 km s _1 , would displace a galaxy by about this distance. 



of a nearby galaxy's superwind: unusually strong metal 
absorption or unusually weak HI absorption. The two 
are contradictory, since metal lines are almost always 
produced by gas with a large HI column density. The 
association of either with superwinds could be strength- 
ened if their presence was correlated with properties of 
the nearby galaxy. 

Figure 1171 compares various properties of galaxies 
with /iMpc < 0.5 (lighter bars) and /iMpc > 0.5 
(darker bars). Star-formation histories and stellar 
masses M* were estimated by fitting model spectral- 
energy distributions to the galaxies' observed photom- 
etry {U n GKJK s +Spitzer-IRAC 3.6, 4.5, 5.8, 8.0/xm for 



galaxies in Q1700; U n GlZJK s for galaxies in Q1623 and 
Q2343; and U n GlZK s for the galaxy in Q2346). The 
models assumed exponentially declining star-formation 
rates with varying time constants and ages, and each was 
subjected to a varying amount of reddening by dust that 
followed a Calzetti (2000) law. See Shapley et al. (2005) 
and Erb et al. (2005, in preparation) for details. Velocity 
differences between the galaxies' interstellar absorption 
lines, Lyman-a emission line, and nebular emission lines 
were calculated as in § [3 The rms widths of the nebular 
lines were estimated by fitting the each line's profile with 
a Gaussian and subtracting the resolution in quadrature. 
Further details can be found in Erb et al. (2004). The 
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Fig. 17. — Characteristics of galaxies with little nearby in- 
tergalactic HI (/iMpc > 0.5). The distributions of these galax- 
ies' properties are shown with dark bars; the distributions for 
the rest of the < Ih^ 1 Mpc NIRSPEC sample are shown in the 
light histogram. The number above each panel shows the approxi- 
mate Kolmogorov-Smirnov probability of observing the data if the 
two distributions were drawn from the same parent population. 
This probability is calculated with an approximation (subroutine 
"probks" in Press et al. 1992) that becomes poor for very small 
sample sizes (as, e.g., in the bottom center panel.) The abscis- 
sae show: r — time constant in the assumed exponential star- 
formation history SFR oc e -t / T ; SFR — star formation rate; M„ — 
stellar mass; E(B — V) — dust reddening, for a Calzetti (2000) dust 
curve; Age — estimated time since beginning of current episode of 
star-formation; Af; a b s — velocity difference between galaxy's near- 
IR nebular redshift and its interstellar absorption lines; At)L y _ a 
— velocity difference between the near-IR nebular redshift and the 
Ly-cr emission line (which was usually not detected); am — ve- 
locity width of the near-IR nebular emission lines. 



data in these figures are available in Table |3J Numbers 
above each panel give the approximate significance of 
the difference in distributions for the two galaxy popu- 
lations, calculated with the routine "probks" of Press et 
al. (1992). Small numbers correspond to more signifi- 
cant differences. None of the differences are especially 
significant given the number of properties examined, al- 
though galaxies with weak nearby HI absorption may 
have somewhat higher star-formation rates. 

Figure HSI is similar, except here the sample is divided 
by the total detected CIV column density within 200 km 
s _1 of the galaxy redshift. Galaxies with Aciv > 10 13 
cm~ 2 are indicated with darker bars, the remainder with 
lighter bars. This differences in this figure may be 
marginally significant. Galaxies near intergalactic CIV 
systems appear to have SEDs that require younger ages, 
shorter time constants, and a wider range of dust red- 
dening. It would be a mistake to interpret this too liter- 
ally: the parameters have large and correlated uncertain- 
ties and our assumed star formation histories are grossly 
oversimplified. The figure suggests, however, that there 
are genuine differences between the SEDs of galaxies with 
and without nearby CIV absorption. The main empiri- 
cal difference is a bluer mean optical-to-near- IR color for 
galaxies near CIV systems, which our SED fitting inter- 
prets as an excess of A stars relative to F and G stars. 



Fig. 18. — Similar to Figure ll7l except here we compare galaxies 
with and without substantial nearby intergalactic CIV absorption. 
Dark bars are for galaxies with Nqiy > 10 cm -2 , where Nqyv ls 
the total detected intergalactic CIV column density (cm -2 ) within 
200 km s _1 of the galaxy redshift, and light bars are for galaxies 
with TVciv < 10 13 cm -2 . 



Missing from figures El an d El is anv indication of 
the large-scale environments that contain the different 
galaxy types. Shot noise prevents us from estimating the 
local galaxy density near any single object, but we can 
estimate the mean galaxy density around an ensemble 
of objects by counting their observed number of galaxy 
neighbors and dividing by the number expected in the 
absence of clustering. Figure El shows the result for the 
different galaxy ensembles discussed above. Galaxies in 
denser environments are more likely to have detectable 
intergalactic CIV absorption within l/i -1 comoving Mpc, 
but local galaxy density does not seem to affect whether 
a galaxy has /im P c > 0.5 or /i Mpc < 0.5. 



7.1. 



7. DISCUSSION 
Superwinds at 2 ^ z ?z 3 

Galaxies would presumably be associated with QSO 
absorption lines even if there were no winds. This subsec- 
tion discusses how strongly (or weakly) our observations 
support the existence of superwinds around the galax- 
ies. We consider below the possibility that intergalactic 
metals were produced at very high redshifts (z ~ 10). 

The blueshifted absorption and redshifted Lyman-a 
emission in high-redshift galaxies' spectra is the strongest 
evidence for winds from the galaxies. It is expected in a 
wide range of models in which the stars are surrounded 
by outflowing gas (e.g., Tenorio-Tagle et al. 1999; Zheng 
& Miralda-Escude 2002) and is not (as far as we know) 
expected in any other class of realistic models. The sta- 
tistical significance is high, since the same pattern has 
been observed in the spectra of dozens of galaxies (e.g., 
figure^. This observation only shows that the outflow- 
ing gas lies somewhere outside the stellar radius, how- 
ever. It does not imply that the outflow will be driven 
into the surrounding IGM. 

The strength of CIV absorption at b — 40 kpc (Fig- 
ure [5J may suggest that the outflows have normally ad- 
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redshift (i.e., fl ~ 1) approximation 
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Fig. 19. — Average large-scale environments for differ ent t ypes 
of galaxies in the < lh' 1 Mpc NIRSPEC sample of figure FT51 The 
environment is quantified through the relative local galaxy density 
p/p, which we calculate by counting the number of neighboring 
galaxies in the full sample within A9 = 200", Av = 500 km s — 1 of 
each galaxy in the subsample and dividing by the number expected 
in the absence of clustering. Galaxies near strong CIV absorbers 
(Nciv > 10 13 cm -2 ) tend to reside in denser environments, pos- 
sibly indicating that the CIV originates outside of the galaxy. In 
contrast, the prevalence of galaxies with little nearby Ly-a absorp- 
tion does not appear to be strongly influenced by environment. 



vanced to this radius, but even so this is only about half 
way to the virial radius (for Mtotai ~ 10 12 Mq in the 
Lambda-CDM cosmology favored by WMAP; Spergel 
et al. 2003). The minimum allowed velocity range of 
the gas at b — 40 kpc (~ 300 km s _1 ) is not large 
enough to guarantee escape. In any case the similarity 
of the CIV absorption strength in Lyman-break Galax- 
ies' (LBG) spectra and at impact parameter b — 40 kpc 
could just be a coincidence; it does not rule out the 
idea that the gas at 40 kpc is falling into or orbiting 
within the galaxy's potential. Although anisotropies in 
the galaxy/ CI V-system correlation function seem to dis- 
favor infall and favor rapid outflows for the origin of the 
CIV, CIV systems tend to be found near galaxies which 
may be too young to have driven winds far into their 
surroundings. The small number of distinct galaxy-CIV 
pairs leaves the situation unclear. 

Measuring the spatial correlation of galaxies and met- 
als on large (Mpc) scales would seem to provide a pow- 
erful way to distinguish between different scenarios for 
intergalactic metal enrichment. If the metals were pro- 
duced in LBGs, they should have a spatial bias 11 similar 
to the galaxies' bias of b ~ 2.5 at z — 3 (Adelberger et 
al. 2005). They would have b — 1 if they were produced 
at any redshift by the numerous dwarf galaxies that are 
lcr fluctuations, and b ~ 1.9 at z = 3 if (as envisioned 
by Madau, Ferrara, & Rees 2001) they were produced 
at z = 9 by galaxies that were 2a fluctuations. These 
estimates of the bias exploit Mo & White's (1996) high- 

11 on scales larger than the maximum wind radius 



for the bias at redshift z of objects that were ^-sigma 
fluctuations at the earlier redshift Zi. The galaxy-metal 
cross-correlation length would therefore be equal to the 
galaxy-galaxy correlation length if the metals were pro- 
duced by LBGs and smaller for the other two cases. Un- 
fortunately the implied bias for enrichment at z ~ 10 
depends sensitively on the unknown value of v. The pa- 
rameter v would not have to be very different from the 
default assumption of Madau, Ferrara, & Rees (2001) to 
make the bias exactly equal to the bias for enrichment 
at z ~ 3. In any case, we have measured the galaxy- 
CIV correlation length, not the galaxy-metal correlation 
length. Our measured cross-correlation function there- 
fore does not have an unambiguous interpretation. 

There are only 9 galaxies in our sample with precise 
(near-IR nebular line) redshifts and little HI within l/i -1 
comoving Mpc, and their statistical significance is even 
lower than one might imagine. There are two reasons. 
First, as discussed above, the measurement is difficult. 
To determine whether a galaxy lies within ^ Ih^ 1 Mpc 
of a narrow absorption line in a QSO spectrum, one needs 
a good understanding of many possible sources of ran- 
dom and systematic error in the estimated redshifts. In 
about half the cases there is no absorption line anywhere 
near the galaxy redshift, but in the remainder the galax- 
ies could conceivably be associated with nearby absorp- 
tion lines if redshift errors were somewhat larger than 
we estimate. Reducing the number of galaxies with little 
nearby HI by 50% would begin to make our observations 
consistent with the predictions of windless SPH simula- 
tions. Second, 5 of the 9 galaxies lie within a single field, 
Q1623, which contains one of the largest known concen- 
trations of QSOs with 2 ;$ z ^ 3. This is a consequence 
of the dense spectroscopic sampling we obtained in the 
field, but it raises the possibility that our result might 
have been different had we surveyed a more representa- 
tive part of the universe. The spatial clustering strength 
of the high-redshift galaxies in this field is certainly not 
typical, for example (Adelberger et al. 2005b). 

Even without its questionable statistical significance, 
the apparent lack of HI near some galaxies would not 
have a straightforward interpretation. This is due to the 
complexity of the interaction between winds and galax- 
ies' inhomogeneous surroundings. Analytic models treat 
this in an extremely crude way, and existing numerical 
simulations are unable to resolve either the shock fronts 
or the instabilities that result when the hot wind flows 
past cooler intergalactic material. As a result it is un- 
clear if winds would destroy the HI near galaxies. Other 
effects that simulations do not resolve (e.g., cooling in- 
stabilities) might reduce the covering fraction of HI near 
galaxies. Even if there were low-density regions near 
galaxies in SPH simulations, they might not be recog- 
nized since the density is generally estimated by smooth- 
ing over the few dozen nearest particles, particles which 
(by definition) are likely to lie outside of any local under- 
density. Finally, photoionization could play some role. 
Kollmeier et al. (2003) and Adelberger et al. (2003) 
have shown that the galaxies' radiation cannot destroy 
all the HI within lh~ l comoving Mpc without overpro- 
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ducing the UV background, but photo-ionizing the weaker 
clouds around one third of the galaxies might be possi- 
ble. It therefore remains unclear if the lack of HI near 
some galaxies is a unique or even an expected signature 
of superwinds. The weak correlations between galaxy 
properties and the strength of nearby intergalactic HI 
absorption may suggest that the HI absorption is unre- 
lated to winds. 

7.2. Metal enrichment at z ~ 10 

Although our observations do not provide unequivo- 
cal evidence for superwinds at 2 ^ z ;$ 3, it would be a 
mistake to conclude that they favor intergalactic metal 
enrichment at very high redshifts. In the few cases 
where the two scenarios have clearly different predictions 
(e.g., anisotropies in the galaxy-CIV correlation function 
on small scales), the observations seem to favor lower- 
rcdshift enrichment. In addition, enrichment at very 
high redshifts seems implausible on other grounds. Fig- 
ure 1201 illustrates a major problem: stars that form be- 
fore z ~ 10 probably produce an inconsequential amount 
of metals compared to the stars that form afterward, 
so their metals could be dominant in the IGM only if 
virtually no metals were able to escape into the IGM 
at later times. In order for 90% of the intergalac- 
tic metals observed at z ~ 3 to have been produced 
at 10 ^ z ^ 15, for example, supernovae ejecta would 
have to be at least 100 times more likely to escape 
their galaxy at z = 10 than at z = 3. This follows 
from the fact that the total intergalactic metal density 
at z ~ 2, p mct ~ Pcr^mct ^ 1.2 x lO 5 M Mpc" 3 for 
Q met ~ 4.4fJ c and fl c ~ 2.3 x 10~ 7 (Schaye et al. 2003), 
would be roughly twice the total metal production at 
10 < z < 15 and five times smaller than the metal pro- 
duction at 2 z ;$ 5 if the cosmic star-formation density 
were roughly constant for 2 < z < 15, the stellar density 
at z ~ 2 were ~5x 10 7 M Q Mpc -3 (e.g., Dickinson et al. 
2003; Rudnick et al. 2003), and 100M Q of star formation 
produced 1M Q of metals. 

Is the metal escape fraction likely to be 100 times 
larger at z ~ 10 than z ~ 3? Two factors are com- 
monly thought to increase the likelihood of escape at 
high redshift: galaxies are less massive and the inter- 
galactic medium lies closer to them than at later times 
when the universe has expanded more. Both are offset by 
competing effects, however, and it is not clear to us that 
very high redshift galaxies have much of an advantage at 
all. 

First consider the galaxy masses. The gravitational po- 
tential energy of a galaxy scales as M 2 /r, or M 5 / 3 , while 
the energy released by a galaxy's supernovae scales as 
fM, where / is the efficiency of star formation. If / were 
independent of M, smaller galaxies would be more eas- 
ily unbound by supernova explosions. This is the stan- 
dard assumption, but in fact / increases with M and 
the situation is not clear cut. For example, compare 
the star-formation efficiencies of Lyman-break galaxies 
at z ~ 3 to those of galaxies at z ~ 10. The ratio of stel- 
lar to baryonic mass for Lyman-break galaxies at z ~ 3 
is M*/M 6 - 0.1 (Adelberger et al. 2004), which implies 
/ ~ 0.1. An upper limit to the star- formation efficiency 
at z = 10 can be derived by assuming that the star- 
formation density is constant for 2 < z < 15. (In fact 
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Fig. 20. — The amount of metals produced in massive galaxies 
at z ~ 2 (LBG) compared to the amount of metals in the IGM at 
z ~ 2 and to the estimated maximum amount of metals produced 
at 10 < z < 15. Star formation at 10 < z < 15 probably does not 
produce enough metals to account for the observed enrichment at 
z ~ 2. Even if it did, the fraction of metals that escape into the 
IGM would have to be ~ 100 times higher at 10 < z < 15 than at 
z ~2-3 for 90% of the intergalactic metals at z ~ 2 to have been 
produced at 10 < z < 15. 



it appears to be significantly lower at z ^ 5; Bunker et 
al. 2004.) In this case the stellar density at z = 10 
will be f2*/f2b = 0.001, scaling from the z ~ 2 stellar 
density estimated by Dickinson et al. (2003) and Rud- 
nick et al. (2003), while the fraction of baryons in the 
halos with v c ^ 10 km s _1 that are able to host star for- 
mation (Dijkstra et al. 2004) is Cl h /Cl b ~ 0.08. That 
implies an upper limit to the star-formation efficiency of 
/ ~ Sl*/£!h ~ 0.01 for typical galaxies at z ~ 10. The 
order-of-magnitude increase in star-formation efficiency 
at lower redshift removes any advantage higher-redshift 
galaxies might gain from their smaller masses. A sim- 
ilar argument shows that the star-formation efficiencies 
of LBGs are far higher than those of smaller galaxies at 
similar redshifts 2 ;$ z ^ 3. 

Now consider the supposed benefit from the small size 
of the universe at z ~ 10. If supernova ejecta at z ~ 10 
extended from their galaxies into the receding Hubble 
flow, their filling fraction at lower redshift would be 
boosted by the subsequent expansion of the universe. 
In fact, however, the ejecta will not normally reach the 
Hubble flow. Instead they are likely to be swept back 
into their galaxy by the ongoing process of structure for- 
mation. To illustrate the point, we show in Figure |^ 
the distance to the nearest halo of mass M ^ IO^Mq 
at z = 2.12 for every particle in the GIF-LCDM sim- 
ulation (Kauffmann et al. 1999) that lay within lh^ 1 
comoving Mpc of a galaxy (i.e., halo) at z = 10 or z = 5. 
The result is stark: the overwhelming majority of these 
particles end up inside the virial radius of a galaxy at 
z ~ 2. Since the stalled ejecta of galaxies at z = 10 
or z — 5 will be swept along by the movements of the 
material that surrounds the galaxies, they should largely 
end up inside galaxies at z ~ 2 as well. This exercise 
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Fig. 21. — Distance to the nearest massive (M > lO n M ) halo 
at 2 ~ 2 for all GIF-LCDM simulation particles whose distance 
r to the nearest halo satisfied 2r v i r < T < 1 comovingMpc at 
z = 10 (top) or z = 5 (bottom). Dashed lines show the particles' 
original (higher redshift) distribution of r/r v i r ; the solid shaded 
histograms show the distribution at z ~ 2. (The nearest halo and 
its virial radius are usually significantly different at the initial and 
final times.) These particles initially lie at larger radii than those 
expected for the metals ejected by very high redshift winds, yet 
they mostly end up inside halos by z ~ 2. This suggests that the 
metals ejected at z ~ 5 and z ~ 10 will generally also lie in massive 
halos at z ~ 2, not in the IGM. 



may be slightly misleading, since the GIF-LCDM sim- 
ulation only resolves halos of mass M £ 10 11 M Q (i.e., 
~ 4cr fluctuations at z ~ 10), not the smaller halos be- 
lieved to be most responsible for polluting the IGM at 
z ~ 10. However, a simple calculation shows that a sig- 
nificant fraction of the metals from the lower-mass pro- 
genitors should also end up inside galaxies at z ~ 2. The 
estimated stalling radius for winds from small galaxies 
at z ~ 10 is ~ 100 comoving kpc (e.g., Madau et al. 
2001), which is the Lagrangian radius for a halo of mass 
1.7 x 10 8 M Q . If these galaxies' typical descendants at 
z ~ 2 have masses significantly larger than this, the met- 
als will likely have been swept inside them. According to 
the extended Press-Schechter formalism, ~ 85% of the 
galaxies at z = 2 that descended from 2a fluctuations at 
z = 10 will have masses that exceed this threshold by 
an order of magnitude (see, e.g., equation 2.16 of Lacey 
& Cole 1993). A substantial fraction of the metals pro- 
duced at 10 < z ;$ 15 should therefore be locked inside 
galaxies by z ~ 2. Once there, some of them are likely to 
cool further, fall towards the center, and disappear from 



view. Since the total metal production at 10 



15 is 



at best comparable to the observed metal content of the 
IGM at z ~ 2 (Figure !2U)l . this suggests that the inter- 
galactic metallicity at z ~ 2 must receive a significant 
contribution from some other source. A corollary is that 
the metal content of the IGM would drain into galaxies 
and decrease over time if it were not continually replen- 
ished. The observed constancy of the IGM metallicity 
(e.g., Schaye et al. 2003) therefore also seems to require 



metals to escape from galaxies at lower redshifts. 12 
8. SUMMARY 

We reviewed the status of our search for direct evidence 
of large-scale outflows around UV-selected star-forming 
galaxies at 2 ;$ z ^ 3. These are our principal observa- 
tions: 

• The gas that lies within 40kpc of LBGs produces ex- 
tremely strong absorption lines (iVciv 3> 10 14 cm -2 ) in 
the spectra of background galaxies and QSOs (Figure[5J. 
The large equivalent widths of the CIV absorption in low- 
resolution spectra imply that the absorbing material has 
range of velocities of at least Av = 260 km s _1 . The 
absorption produced by this gas is similar to the inter- 
stellar absorption seen in LBGs' spectra, suggesting the 
LBGs' outflowing "interstellar" gas may actually lie at 
radii approaching 40kpc. 

• For roughly half of the LBGs, CIV absorption lines 
with column density New ~ 10 14 cm~ 2 are observed 
out to impact parameters of ~ 80 kpc (Figures |S] and 
This implies that roughly one-third of all "intergalactic" 
absorption lines with iVciv ~ 10 14 cm -2 are produced by 
gas that lies within ~ 80 kpc of an LBG. Galaxies too 
faint to satisfy our selection criteria could easily account 
for the remainder. In some cases the absorbing gas has 
substructure on half- kpc scales (Figure 0- 

• The cross-correlation function of galaxies and CIV 
systems with New ~ 10 12 5 cm" 2 appears to be the same 
as the correlation function of galaxies (Figure ITT)). This 
implies that CIV systems and galaxies reside in similar 
parts of the universe and is consistent with the idea that 
they are largely the same objects. We also find a strong 
association of OVI systems with galaxies (Figurc lToT) . On 
small scales the redshift-space cross-correlation function 
is highly anisotropic (Figure lT^|l . The metal-enriched gas 
must therefore have large velocities relative to nearby 
galaxies. The required velocities appear to exceed the 
galaxies' velocity dispersions but are similar to the galax- 
ies' observed outflow speeds. 

• In contradiction to the earlier result of Adelberger et 
al. (2003), we find that the gas within lh' 1 comoving 
Mpc of LBGs usually produces strong Ly-a absorption 
in the spectra of background galaxies (Figure I14fl . The 
absorption is weak (mean transmitted flux within l/i -1 
Mpc of /iMpc > 0.5) in only about one case out of three 
(Figure IT5|l . Even so, the weakness of the absorption in 
these cases remains difficult to understand. Since high- 
redshift galaxies reside in dense parts of the universe, 
with large amounts of hydrogen and short recombina- 
tion times, one would expect them to be surrounded by 
large amounts of HI. The HI has presumably collapsed 
into clouds and sightlines to the background QSOs might 
occasionally miss every cloud near a galaxy, but the SPH 

12 If one removes the assumption that the metals injected into 
the IGM at a given redshift must lie near galaxies that existed 
at that redshift, then it is possible to find a spatial distribution 
of intergalactic metals that matches intergalactic absorption-line 
statistics and does not drain into galaxies at a significant rate for 
2 z ;$ 4 (e.g., Schaye et al. 2003). The removal of this assump- 
tion seems dubious to us, however. Metals could exist in parts of 
the IGM that were never near galaxies only if they were produced 
by Population III stars, presumably, and the amount of Population 
III star formation would have to be inordinately large t o m atch the 
observed intergalactic metallicity at z ~ 2 (e.g., FigureCT see also 
Aguirre et al. 2004). 
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simulations of Kollmcicr et al. (in preparation) suggest 
that the chance of this is very low (Figure EJ- 

• We were unable to identify any statistically signif- 
icant differences in age, dust reddening, stellar mass, 
kinematics, or environment between galaxies with weak 
nearby HI absorption and the rest, although galaxies 
with weak absorption may have higher star-formation 
rates (Figure ITT|) . Galaxies near intergalactic CIV sys- 
tems appear to reside in relatively dense environments 
(Figure EJ) and to have distinctive SEDs characterized 
by blue colors and young ages (Figure ITB)) . 

As discussed in § none of these observations provide 
unequivocal evidence for superwinds at z ~ 2-3. This 
is mostly because it is unclear how superwinds would 
affect the correlation between galaxies and intergalac- 
tic absorption lines. Metals and galaxies would lie near 
each other at z ~ 3 even if there were no active outflows, 
for example. We argued in § that star-formation at 
10 < z < 15 may not produce enough metals to account 
for the intergalactic mctallicity at z ~ 2-3, and that 
many of these metals would be buried inside galaxies by 
z ~ 2-3, but our simple arguments need to be checked 
with numerical simulations. The interpretation of our 
data will remain ambiguous until that time. The argu- 
ments presented in § should make it clear that metal 
enrichment at 5 £ z £ 15 will be much harder to rule out 
than metal enrichment at 10 z ^ 15. 

There is only one case where the data seem easy to 
interpret and we are limited by the small size of our 
sample. This is the redshift-space distribution of met- 
als within ~ 200 comoving kpc of galaxies. Mapping the 
velocity offsets between galaxies' stars and their metals 
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as a function of impact parameter would help show if the 
detected metals are flowing out of, falling into, or orbit- 
ing within the galaxies' halos. The velocity offsets can be 
mapped by obtaining higher resolution spectra of more 
close pairs (e.g., the right panel of figure or by obtain- 
ing more near-IR nebular redshifts for galaxies near the 
QSO sightline (e.g., Figure Et - That seems the sensible 
way to proceed. 
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